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Abstract 

New spectral representations for f n and chiral condensate are de- 
rived in QCD and used for calculations in the large N c limit. Both 
quantities are expressed in this limit through string tension a and 
gluon correlation length T g without fitting parameters. As a result one 
obtains (qq) = —N c a 2 T g ai, f n = \pST c oT g a<i-, with a\ = 0.0823, 02 = 
0.30. Taking a = 0.18GeV 2 and T g = 1 GeV , as known from ana- 
lytic and lattice calculations, this yields (qq)(n = 2GeV) = — (0.225GeV) 3 , 
f n = 0.094 GeV, which is close to the standard values. 

1 Introduction 

The Chiral Symmetry Breaking (CSB) is known to occur in QCD at large 
N c , if confinement is preserved in this limit pQ. Lattice calculations for 
N c — 2, 3 indicate that confinement and CSB coexist in the confinement 
phase at T < T c and disappear simultaneously above T c At larger A^ c it 
was found on the lattice that the 1/N C corrections to all observables studied 
are not large [H], suggesting that a smooth limit at large A^ c is possible. 

In the framework of the Field Correlator Method FCM [1] the dynamics 
of confinement and deconfinement is associated with the set of field correla- 
tors D ( i ^ vl 4XnVn (x 1 ,...x n ) = (F^^x^.-.F^^Xn)) 1 of which the lowest one 

1 parallel transporters are here omitted for simplicity 



D^ 2 \xi,X2) = D^tyXi — X2) plays the dominant role Moreover, D( 2 \x) 
was calculated on the lattice jH] and its confining part, D(x), was shown to 
disappear exactly above T c 

In jSJE] also CSB was found as a consequence of confinement and in PITU] 
the Effective Chiral Lagrangian (ECL) was derived from the 4q interaction 
term using D( kernel. 

The resulting ECL in [TU] has a general structure which can be reduced 
to the expressions derived in the framework of the instanton model ^1] or 
the NJL model when the corresponding kernels are introduced there. 

In the case of confinement, the effective quark mass operator M(x) in 
QCD obtained in [§J QUI contains the effect of the scalar confining string 
connecting the quark to the nearest antiquark. Moreover all invariant quark 
Green's functions can be expressed at large N c through the string spectrum 
as it was done in in the PS channel. 

The phenomenon of CSB was shown in EH] as occurring due to the 
spontaneous creation of the scalar string (similar to the creation of the scalar 
condensate in nonconfining models jT^. ) which generates CSB and chiral 
Nambu-Goldstone (NG) fields (see eqs. (50-54) in and eqs. (21-24) in 

Since confinement is present in our formalism (in the form of M(x)) one 
can ask the question how confinement fits in the chiral picture of NG spec- 
trum, and in particular how CSB modifies the lowest PS states computed 
in FCM (or in any quark model) taking into account confinement and disre- 
garding CSB. Two such lowest states, and its first radial excitation ir^ 
with masses m(ir^ ') = mo = 0.4 GeV and m(ir^) = mi = 1.35 GeV have 
been computed in FCM, see Appendix 2 below in this paper. It was shown in 
[TUj that the ECL obtained there with account of confinement, has a remark- 
able property: the PS spectrum of confinement transforms due to CSB in 
such a way that vr^ -* becomes a NG pion with the mass satisfying Gell-Mann- 
Oakes-Renner (GOR) relation while the first radial excitation shifts only 
slightly. 

In deriving that property it was essential that all basic quantities in the 
ECL and in particular the pion self-energy operator can be expressed as 
a spectral decomposition in the confinement (string-like) spectrum states, 
which is possible in the large N c limit. 

In this paper we follow this line to obtain a more fundamental relation, 
namely, to calculate the quark condensate (qq) and the pion decay constant 
f n using new spectral representations for these quantities. Since in the latter 
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all masses and coupling constants are expressed via D(x), i.e. via the string 
tension a and the gluon correlation length T g , we have an expression for 
(qq) and f„ in terms of a and T g . The most important role in the spectral 
representations is played by the lowest PS meson 7r(°) - the "to be pion" - 
which is the quark model analog of the pion with mass mo shifted by the 
hyperfme interaction from the p-meson mass. In Appendix 2 we derive the 
mass mo and the corresponding wave-function in the framework of FCM in 
terms of a and a s . 

Having established the connection of (qq) and f n with a, T g stated in 
abstract above, and explained in the text below, it is easy to understand 
that at the the deconfmement transition when a vanishes at T = T c , also 
(qq) and vanish in agreement with lattice data j2]. 

Some specification with respect to the notion of "magnetic confinement" 
[Hj is needed at this point since magnetic counterpart of D(x) and the cor- 
responding spacial string tension stay nonzero above T c . This topic will be 
studied elsewhere. 

The paper is organized as follows. In the next section the ECL is written 
down together with the appropriate expressions for (qq) and f n . In section 3 
the spectral representations for these quantities are derived, with coefficients 
depending on eigenfunctions of the qq system in the pseudoscalar channel. 
Section 4 is devoted to the discussion of results in comparison to lattice 
data and to the concluding remarks. Four appendices are included in the 
paper, containing respectively the evaluation of M(0), derivation of spectral 
representation, explicit calculation of eigenvalues and eigenfunctions in the 
pseudoscalar spectrum, and the contribution of the small-distance region. 

2 The Effective Chiral Lagrangian 

The quadratic part of the ECL for pions was derived in and has the form 

N r - <7( 4 )t 

^ (2) (0) = yJ Mk)M-k)N(k)j^y 4 (1) 

where notations of [10J have been used, a = y 2 - and 

N(k) = \[G^ MM \k) + tr(AMs)} = (rr£ + k 2 )^- + Q(k% (2) 



3 



G (MM ){k) = _ tr{X{y,x) lb M s {x)K{x,y) lb M s {y)y k ^d^\x - y), (3) 



A(x,y) = @ + m + M s )-i. (4) 

As it was shown in JU], two terms in the square brackets in (J2J) cancel 
for k 2 = m = and one obtains the GOR relation for the pion mass J3] 

1 _ 

mN c trA = m\(ipip) M \ = -{m u + m d )\(uu + dd) \ = mlf 2 . (5) 

To calculate the quark condensate, defined in the Minkowskian space 
time, one can write (ipip)M — —N c trA, and use identical transformation 

trA xx = tr( —(M s + m- d) — ) = 

{M s + m + d) '{Ms + m-d) 

(tr( l5 A(x,y) l5 (M s +m)A(y,x)))d*y = -J G^ M \x,y)d 4 y = -G^\k = 0). 

(6) 

Hence trA xx reduces to the zero-momentum component of the qq Green's 
function in the PS channel, which differs from (jSJ) only by vertex operators. 

To define f^, one needs the first term in the k 2 expansion of G^ MM \k) 
Eq.flSJ), (c.f. Eq.Q) so that one has 

G ( MM \k) - G (MM) (0) = + 0(k 4 ) (7) 

2,N c 

As it was argued in [10], both G^ MM ^(fc) and G <yM \k) have spectral represen- 
tations in the large N c limit, with the same set of poles m n , n — 0, 1, 2, m = 

m(7r(°)), 

oo (JM)\2 oo (A/) 

G (MM) (fc) = _ E JC^ = ( 8 ) 

^0 fc + m n ^0 k2 + m l 

In the next section we shall determine the coefficients c n , c^f 1 ' and m„ for 
the lowest states, and now we define (qq) and /„. in terms of spectral sums 
(JHJ). From © and (J7J) one has 

oo (M) oo ( JM)\2 

<^> M = -2iV c £^, /. 2 = 2iV c E^- 0) 

n=0 " l n n=0 '"n 
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The coefficients c n and differ by the presence of the vertex operator 
Ms = M(0) in the latter which is a constant computed in Appendix 1, 
therefore one has = M(0)c n , and limiting oneself to the first term in the 
sum © one obtains 



K^)M|>m^-^ = ^£. (10) 



Ch 



M(0) 



Inserting M(0) = 148 MeV from Appendix 1 and |^| = (225 MeV) 3 , 
f n = 94 MeV, one obtains mo = 437 MeV which is close to the value mo = 
400 MeV calculated in Appendix 3. On the other hand, the sum @ for {ipip) 
is converging more slowly than that for f n , and therefore one has inequality 
in (fTUj) due to the presence of higher terms in (ipip) . 



3 Calculation of (qq) and f n 

The integration region in the space-time integrals in Q, (0 can be split 
in two parts: \x — y\ > T g and \x — y\ < T g . In the first (long distance) 
region the relativistic local potential-type dynamics sets in at space-time 
distances exceeding T g [TH], [IE] and the result can be expressed in terms of 
the spectrum, as will be done below in this section. The second region can 
be treated in the OPE formalism J7j and is considered in the Appendix 4. 
It is shown there that the contribution of this region is parametrically small 
in the parameter aTg -C 1. Only the long-distance contribution is calculated 
in this section below. 

We start with the calculation of and to this end we write the qq Green's 
function G^ MM \k) in terms of cm. and relative coordinates as follows (an- 
other derivation is given in Appendix 2): 

-G^ IM \k) = J d 4 XG (MM \r 12 = 0; R = 0, r' 12 = 0, R' = X, T)e* kX = 

= M 2 (0) V \cp n (0)\ 2 ^dTd 3 ^e- E ^ T - lPX+lkX = 
(27iT 

= M 2 (o)$>n(o)| 2 r e -^ k ^r = M 2 (o)E (11) 

n Jo n ^jm 2 n + k 
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Expanding (jlljl in k 2 and comparing to (J7J) one finds / 2 , 

12 



„_(] m n 



fl = N c M\0) £ (12) 
Comparing Q and (|12j) one finds 



ci M) = ^M(0)^(0), c n = ^>„(0). (13) 
In a similar way one computes (qq) from (jHJ and finds 

- (qq) = N c M(0) £ (14) 

n=0 m « 

Here ip n (r) is the 3d spin-singlet wave- function of qq system, as obtained 
e.g. in the relativistic Hamiltonian method of FCM ^3], or else in the Bethe- 
Salpeter equation with the kernel not depending on relative time, as it is 
discussed in Appendix 2. 

The accuracy of the method with respect to calculation of |<£> n (0)| 2 can be 
checked by comparison of predicted leptonic width with experiment, as it is 
done in Appendix 3. Taking into account both color Coulomb and confining 
interaction one has 

, (0)l , = ^+pm (15) 

where \i n is the constituent energy (mass) computed through a [TS]; we refer 
the reader to the Appendix 3 for the details of calculation of m n and |y? n (0)| 2 . 

As it is shown in [TH] and discussed in Appendix 2 , the masses m 2 ,/z 2 
grow linerly with n in the large iV c limit, hence the sum ()14j) for (qq) is 
formally diverging if the spectrum of radially excited mesons extends to in- 
finitely large masses. In fact the experimental spectrum can be followed up 
to the mass values around m cont = 2.5 GeV, where resonances become very 
wide and strongly mix between themselves and with hybrids, forming the 
continuum of states. Following the ideology of the QCD sum rules [T7j one 
could replace this continuum by the perturbative diagrams, which do not 
contribute to (qq). Therefore we shall keep the first 3 terms in the sum ()14|) 
over n (the term with n = 3 gives negligible contribution). As was men- 
tioned in the begining of this section the relativistic potential description of 
G( MM '(ri2,R,0,ri2,R',T) is possible only for the time T > T , while for 
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T < T , T ~ T g one should use the properties of the qq Green's function 
G, as given by the OPE [T7j. As it is discussed in Appendix 4, the region 
of small times and relative distances covered by the OPE treatment, gives a 
contribution to (qq) proportional to 0(a 5 ' 2 7&, am) and therefore can be dis- 

9 

regarded for light quarks and small T < T g . As a result one should exclude 
from the integration over dT in (jllj) the region (0, T ) which brings about 
the following factor in ()14|) instead of l/m n 

I e -m n T 



m n m r . 



(16) 



and in (fHJ) 



I e -m n T 

mi mi 



l + m n T ). (17) 



Keeping for (qq) the first 3 terms in the sum (JUjl and 2 terms in (|12|) one 

has 

{ mg mi m 2 j 

fl = N c M(0) { ^^^ (l + m T ) + ^ (0)e 3 " mir ° (l + mi T )\ . (19) 
I m'Q ml J 

Using (|S!251> , (BOj) one has 

2/ , 0.109GeV 3 2/ „, 0.097GeV 3 2/ „, 0.115GeV 3 
Vo(0) = ^^ , ^i(°) = ^^ . ^ (0) = ^^ • (20) 

m = 0.4GeV, m t = 1.35GeV, m 2 = 1.85GeV. 

For a reasonable estimate we put T = T g = lGeV -1 and the value 
M(0) = 0.148 GeV from Appendix 1, and obtain. 

- (qq) = (0.195GeV) 3 , f w = 0.094GeV. (21) 

One can check that the behaviour of ()18j) for (qq) at small To is smooth, 
e.g. when changing T from T g = 1 GeV -1 to 0.5 T g , the result changes by 
roughly 10%. 

To check the sensitivity to the change of T g , we have taken T g = 1/1.5 
GeV -1 and recalculated all quantities, e.g. from (jA.7J) one has M(0) = 0.12 
GeV. The resulting values are not much changed from (|21j) . 

- (qq)(T g = -^-GeV" 1 ) = (0.189GeV) 3 , f n = 0.076GeV. (22) 

J. .0 
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It is remarkable that f w in (|21j) is very close to the value obtained from the 
pion decay and used in the chiral perturbation theory [TH] f n = 93 MeV. At 
the same time \(qq)\ is somewhat less than the standard value (240 MeV) 3 , 
and we discuss in the concluding section section the scale dependence and 
comparison to existing lattice measurements. 

4 Discussion and conclusions 

The quark condensate and f n are given by Eqs.(fT%|). (|T9*|) and (J22J), where all 
quantities can be expresses through m, a and T g , since (Pn(0), m n and M(0) 
are expressed through these quantities, while T can be taken in the region 
of plateau and e.g. equal to T g . In this way one obtains (m = 0, o = 0.18 
GeV 2 , T = T g = 1 GeV-\ and Eq.(j22D for (qq)) 

f n £ 0.094GeV, - (qq) £ (0.20GeV) 3 . (23) 

Several corrections should be added to this results. First of all, the short 
distance contribution to (qq) is of relative order y/aT g ~ 0.45 and can sub- 
stantially increase the result. Another essential point is the value of T g , 
which increases in the presence of dynamical quarks, and can be smaller if 
gluelump data (TOJ are instead taken into account, T g = 0.7 GeV -1 . This 
influences significantly the value of M (0), however an independent check can 
be made since M(0) also enters the strong decay matrix element, and the 
value M(0) = 0.148 GeV is reasonably close to the phenomenological value 
known from the 3 Pq model [20J. 

We are now in position to compare (|2"3*|) with the lattice data. There the 
computation was done in the quenched case for Wilson fermions [21] and also 
for the overlap action [22J. Before using the evaluation coefficient for (qq), 
one can compare the result ()23|) which does not contain any scale /i, and any 
evolution corrections, with the so-called Renormalization Group Invariant 
(RGI) lattice measurements, which yield j2Ij 

- (qq)™ 1 = [(206 ± 44 ± 8 ± ±5) MeV] 3 . (24) 

This value is in reasonable agreement with (J23|) . As the next step we take 
the evolution coefficient for (qq) computed in [22- ip>f — 0, N c — 3) 

CfS(//) = [a s (/i)]- 4/11 {l - 0.219a s - 0.1054a s 2 }. (25) 
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For \i = 2 GeV taking a s ~ 0.3, and identifying (qq) in ()23|) with {qq) RGI 
one obtains for the long-distance contribution to the condensate 

(qq)(n = 2GeV) S (qq) RGI Cj^ = -(225MeV) 3 . (26) 

This value, given in the abstract of the paper, is obtained without inclu- 
sion of the coefficient used on the lattice j2Ij to calculate the transition from 
the lattice RGI result to the MS scheme, this coefficient is anyhow close to 
unity. 

The lattice value at \i = 2 GeV for Wilson quarks in |2T)] 

(qqy^fa = 2GeV ) = -K 242 ± 9)MeV)] 3 (27) 

and differs from the result [22]: - (282 (6) MeV) 3 ( 176 a 6 "^ eV ) 3 . An independent 
estimate from the QCD sum rules yields [24 

(qq)(ti = M N ) = -[(225 ± 9)MeV] 3 . (28) 

As a result one can see that our long-distance contribution to (qq), Eq. (|2EJ), 
is somewhat smaller than the lattice data ()27j1 . but is certainly in the same 
ballpark, and the evaluation of the short- distance contribution is important 
to improve the accuracy of calculation. 

At the same time the resulting value f m (}23|) is in good agreement with 
the standard value, obtained from the pion decay and used in the chiral 
perturbation theory [THj . 

The method used above can be easily applied to the case of nonzero quark 
mass m and the SU(3) flavour group to calculate (ss), fx etc., which will be 
published elsewhere [2*5] . 

The financial support of INTAS grants 00-110 and 00-366 is gratefully 
acknowledged. 
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Appendix 1 

Calculation of the vertex mass M(0) 

One starts with the definition of the nonlocal mass operator Ms(u,v), 
given in PCU1 (see e.g. Eq.(24) in PU) 

M s (u, v) = (7 M A(w, v)-f^) sc J(u, v). (A.l) 

The mass operator enters in the gauge-invariant Green's functions, see 
e.g. Eq.(3), via the quark propagator A(x, y), Eq. (4), where M s (z, u) enters 
at all intermediate points, and also at initial and final points and ^^i 
where the nonlocal pion x) is emitted. According to the prescription 
given in [10J, we choose the set of contours C(z) for all intermediate points z 
in the Green's function C7(^7p, ^^), which minimizes the mass eigenvalues. 
One simple choice is to take the contours C(z) from z a long the shortest 
way to the x\ axis passing through and and along 24 axis to the 
origin at the point ^rp. 

When Ms is situated at the initial or final point of the qq Green's function, 
i.e. at the points M s (x,x) or M s (y,y), where the qq or pion is created or 
annihilated, then it is convenient to choose points x, x on the axis 1 with the 
origin at ^rp. In this way one obtains for x\ > 0, X\ = y\ < 0, y^ = £4. 

J(x,y) — / du I dvD(u — v, x 4 — y 4 ). (A. 2) 

JO Jy x 

It is convenient to use for D the Gaussian form, 



£>(x,x 4 ) = £K0)e W 

which yields 

J(x,y) = -e (1-e ~^~). (A.4) 

7T 

Now one has to estimate the scalar part of the quark Green's function 
A(x, y) in (|A.1|) . for which in jHj it was found that it behaves as a smeared 8- 
function with the smearing radius equal to 1/a/o" (see Eq.(24) in the second 
ref. of 8j). We simplify this form taking 

/ fT \3/2 r 

A{x, y) = {-) *■ j A{x, y)rf 3 (x - y) = 1. (A.5) 



a 



(A.3) 
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To obtain the localized form of the vertex function 



M s (x)=M(p) = I M ( i\x 1 y)d\x-y) 



(A.6) 



we substitute in (IA.6J1 J(x, y) from (jA.4j) and A(x, y) from (|A.5|) to get finally. 



M(0) 



' a 4 T 2 ' 



3/2' 



2<tT„ 



2ctT„ 



7T 



(A.7) 



one obtains M(0) 



i.e. exactly the value 



In the limit o~T 2 

which appears in the strong decay vertex of the string in [20] ■ This is not 
surprising since in both cases M(0) is a mass corresponding to a piece of the 
string with the length of the order of T g , hence M(0) ~ crT 9 . The factor rj 
in (jA.7|) describes the attenuation due to the nonlocality of A(x, y)) at small 
|x — y| for light quarks. For heavy quarks this factor tends to zero since 
the localization of A(x,y) becomes more strong, indeed the quark Green's 
function A for m — > oo is proportional to 6^ (x — y ) see |S] . Effectively for 
nonzero m this can be described replacing 77 in (jA.7|) by the factor 



i] — > 7](m) 



{oT g + m)AT g 
(aT g + m)4T g + 1 



3/2' 



(A.8) 



For light quarks and o = 0.18 GeV 2 , T g — 1 GeV 1 , the factor r\ = r/(0) 
is 1 - 0.27 = 0.73, and from (1A~71) one gets M(0) w 0.148 GeV. 

Appendix 2 

Derivation of the spectral representations, Eqs. (TT2I) and ( IT31) 

Consider the gg Green's function of the type given in Eq. (jUJ) 

G r (x, y) = (trrA(x, y)TA{y, x)) (A.9) 

where T = 75, 7^, and the qq Green's function in the 4x4 spinor repre- 
sentation, 

G {qq ) al3r/ 8(x, x; y, y) = (A a/3 (x, y)A jS (y, x)). (A.10) 

Following the standard procedure from one can introduce the cm. 
and relative coordinates, e.g. 



X 



x + x 



Y 



y + y 



x — x } r = y — y 



(A.ll) 



11 



and define 

G^\x,x;y,y) = J d'Pe^ x - Y ^e-^G P (r, r', e, r' ) (A.12) 

Gp satisfies an equation 

(E~E 2 -H 1 )(E 2 -H 2 )G P = p 1 p 2 i (A.13) 

where E = E x + E 2 = P , E 1 - E 2 = 2e, and 

Hi = rriiPi + pcti + (3iM S - (A.14) 

At this point one can exploit the property of Hi that it does not depend 
on relative time tq, and therefore one can integrate in (jA.12|) over de with 
the result HH |27] 



Gj 



oc 



, _ f -P 1 p 2 l — =i p 1 i -l 1 H = Hy + H 2 . 

r o-U J-oo {E x - H 1 )(E 2 - H 2 ) E-H 



r' = 



(A.15) 

As the result one obtains 



fd\X-Y)G^\x,x;y } y) ro =r>=o = (r ,, , 

(A.16) 

One can now express Gp, with T = 75, 



/^(x-y )Gr (x,,) = i £«5MM (A . 17) 

n =0 

where we have defined the relativistic wave- function ip n (r) ~ 7 5 (r|n), ip^ ( r ) ~ 
/?i/?2 (tz I r) , satisfying the Hamiltonian equation 

FVn(r) = £ n Vn(r). (A. 18) 

As it is known from dynamical calculations with the Bethe-Salpeter equation 
[2*H] with the scalar confining kernel the dominant role in ip n (r) is played by 
the 1 So component y2 n (r) of the wave-function, which satisfies the relativistic 
Schroedinger equation with the hyperfine interaction, discussed in the Ap- 
pendix 3. Therefore one can identify ip n (r) — > <p n (r),E n — > m n and the Eq. 
© with the help of (IA~TY|) goes over into Eq. (fTijl . 
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Appendix 3 

Calculation of the masses m n and (p n (0) through a in FCM 

~ p) . (tt) 

The mass eigenvalue m n of the spin-averaged state m " for L = 
can be written as [121 HE] 

m n = M (n,0) + A 5 £ + A c , (A.19) 

where M (n, 0) is the eigenvalue of the spinless Salpeter equation, which can 
be written as 



M (n, 0) = 4/i (n) = 4<^p 2 + m 2 ) n0 . (A.20) 

For m = 0, /io(n, 0) is expressed through -^/a and dimensionless coeffi- 
cients a(n) - zeros of Airy functions [Tl 

3/4 



an 



^ ( n ) = ^ f _i_>J , a (0) = 2.338, a(l) = 4.088. (A.21) 

Taking into account nonzero m one finds Ho(n) from the equation 

TO 2 cr 2/ ^ 

l = _ + - 1 ^a(n). (A.22) 



For large m 3> \At the solution of (|A.22|) is 



m 



1 + 



a(n) ( 2cr fj, \ 2 ^ 3 



(A.23) 



3 \ v m 2 (m + (j) J 

The term A$e is the self-energy correction ^H] which can be written as 

. , Aari(m) 

and 77(777.) is computed through m, for m = it is ?7(m = 0) = 0.9 -j- 1. 

Taking all contributions into account one obtains for the light quarks 
(m = 0) 

rn = 0.652GeV, to 2 , = rh 2 + tt n, n = 0, 1, 2 .. (A.25) 

where Q is computed solely through a and is equal VLq — 1-6 GeV 2 , which 
is close to the experimental slope Q exp (L = 0) = 1.64 ± 0.11 GeV 2 , see jH 
for refs. and discussion. 
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Now we take into account the hyperfine interaction which produces the 
HF splitting between p and 7r states. 



^ F = A HF + A HF , A HF = 2 , r • A.26 



Here i? n (0) = V47T(y9 n (0) is the radial meson w.f., which can be found also 
from the leptonic width of p meson. One has 

iTWnM2 4 / -2\n f 0.091GeV 3 , a s = 

|^o(0)| 2 = ^ (0)(a + -«.<r 2 » = . 109GeV3 ; a = . (A.27) 



These values can be checked vs the leptonic width of p, T e 
while T e * p e - = (6.85 ± 0.11)KeV. 



6.36KeV, a s = 
7.62KeV,a, = 0.3 



Thus one obtains Agf* from (|A~26j) . Af ert 



0.26GeV,a s = 
1HF ~ \ 0.24GeV,a s = 0.3. 
The nonperturbative part A^ F is expressed through the correlator D(x) 
[2Hj and depends on the accepted value of G 2 = ^(F^F^J) , 



We take two values of G 2 = Gff = 0.012 GeV 4 and G 2 = 2G s 2 t . Thus one 
obtains for the lowest mass of PS state in the qq approach (no chiral effects) 

m = m -^A HF = [0.652- ^(0.3^0.35)]GeV = (0.39-0.43)GeV. (A.29) 

As a result we accept the following values for m and m\ (the latter is 
calculated in the same way using (jA.25|) and A HF {n = 1)) 

m = 0.4GeV, m x = 1.35GeV. (A.30) 

Appendix 4 

Small distance contribution to (qq) 

To separate the small-distance contribution we start from Eq. © where 
we take into account the nonlocal structure of Mg(u,v) and put m — 0, 



trA xx = I d%d A ytr( l5 A(x, yh 5 M s (y, u)A(u, x)). (A.31) 
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In the limit when one keeps the most singular part - the free part of A(x, y), 
one has 

A(x, y) « A /r «(x, y) = — 2 Ij-JL; + J + ... (A.32) 

From (jA.lj) one can derive the behaviour of Ms(y,u) at small |y — u| < 

To, T g , 

M s (y, u) ~ ^c|y - w| 2 A(x, y), (A.33) 

where the coefficient c is of the order of unity. 

The nonperturbative part of A(x,y) is not singular (modulo logarithms) 
and proportional to a 3 / 2 , c.f. Eq.(A.5) (apart from the OPE part of A which 
has m(qq) and (F 2 ) terms and even less singular at small x). Finally inserting 
flA~33l . (TP2j) into (lA~3T| and integrating in the region \x— y\, \y— u\, \u—x\ < 
T one can write the short- distance contribution to (jA.31|) as 

trA 3;:c (small distance < T g ) = 0{am, {(r 5/2 -^). (A.34) 

9 

As it is seen from fllHj) the long distance part is 0(a 2 T g ) and is dominant at 
aT 2 ^0,T <Tg. 
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